arXiv:1509.00492v2 [gr-qc] 5 Dec 2015 


Cosmological constraints on superconducting dark energy models 


Zoltan Keresztes^’^,* Laszlo A. Gergely^’^,'^ Tiberiu Harko^,^ and Shi-Dong Liang"^^ 
^Department of Theoretical Physics, University of Szeged, Tisza Lajos krt 84-86, Szeged 6720, Hungary, 
‘^Department of Experimental Physics, University of Szeged, Dom Ter 9, Szeged 6720, Hungary, 
^Department of Mathematics, University College London, 

Gower Street, London, WCIE 6BT, United Kingdom, and 
State Key Laboratory of Optoelectronic Material and Technology, 
and Guangdong Province Key Laboratory of Display Material and Technology, School of Physics and Engineering, 
Sun Yat-Sen University, Guangzhou 510275, People’s Republic of China 


We consider cosmological tests of a scalar-vector-tensor gravitational model, in which the dark 
energy is included in the total action through a gauge invariant, electromagnetic type contribution. 
The ground state of dark energy, corresponding to a constant potential G is a Bose-Einstein type 
condensate with spontaneously broken U(l) symmetry. In another words dark energy appears as a 
massive vector field emerging from a superposition of a massless vector and a scalar field, the latter 
corresponding to the Goldstone boson. Two particular cosmological models, corresponding to pure 
electric and pure magnetic type potentials, respectively are confronted with Type lA Supernovae and 
Hubble parameter data. In the electric case good fit is obtained along a narrow inclined stripe in the 
Qrn — Llv parameter plane, which includes the ACDM limit as the best fit. The other points on this 
admissible region represent superconducting dark energy as a sum of a cosmological constant and 
a time-evolving contribution. In the magnetic case the cosmological test selects either i) parameter 
ranges of the superconducting dark energy allowing for the standard baryonic plus dark matter or 
ii) a unified superconducting dark matter and dark energy model, additionally including only the 
baryonic sector. 

PACS numbers: 98.80.-k, 98.80.Jk, 98.80.Es, 95.36.+x 


I. INTRODUCTION 

The ACDM (ACold Dark Matter) cosmological model, 
in which our Universe is dominated by two major com¬ 
ponents, dark energy and dark matter, has successfully 
explained most of the recent cosmological observations. 
The Planck satellite data from the 2.7 degree Cosmic Mi¬ 
crowave Background full sky survey [1, 2] provided an¬ 
other important recent confirmation. The presence of a 
simple cosmological constant A in the Einstein gravita¬ 
tional field equations provides the late-time acceleration 
of the Universe [3], a phenomenon representing one of 
the most important challenges for present day science. 
The inclusion of the cosmological constant in the gravi¬ 
tational field equations however is highly unsatisfactory, 
since it lacks a well established physical meaning or inter¬ 
pretation, and in addition its numerical value cannot be 
predicted from present day physical knowledge. From a 
theoretical point of view it is desirable to model the mys¬ 
terious dark energy component of the Universe, which 
determines the acceleration of the distant supernovae, 
either as a specific physical field, or as a modification of 
the gravitational force itself. 

An intriguing possibility in this framework is that at 
large, cosmological scales the very nature of the gravi¬ 
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tational interaction changes, such that general relativity 
cannot describe at these scales the dynamical evolution 
of the Universe. Several theoretical models, based on 
various modifications of general relativity, like f{R) grav¬ 
ity (with R the Ricci scalar) [4], f{R,Lm) type models 
(with Lm the matter Lagrangian) [5], the f{R,T) modi¬ 
fied gravity type models (with T the trace of the energy- 
momentum tensor) [6], the Weyl-Cartan-Weitzenbdck 
(WCW) gravity [7], hybrid metric-Palatini f{R,TL) mod¬ 
els (77. being a Ricci scalar formed from a connection 
independent of the metric) [8], f{R,T,R^,j,T^i, gravity 
(with R^ii, the Ricci tensor and the matter energy- 
momentum tensor) [9], or the Eddington-inspired Born- 
Infeld theory [10] have been extensively investigated. For 
a recent review of the generalized curvature-matter grav¬ 
itational models of / {R, Lm)) and f{R, T) type see Ref. 
[11]. Modified gravity models can generally explain the 
late acceleration of the Universe, without the need of 
dark energy, and can also offer some clues on the nature 
of dark matter. 


In a field theoretical approach, among the most impor¬ 
tant dark energy candidates are scalar fields, which nat¬ 
urally arise in particle physics and string theory. Dark 
energy models involving time-dependent scalar fields pro¬ 
vide a cosmological dynamics in accordance with obser¬ 
vations. Moreover, the early evolution of the Universe 
in its inflationary era can also be understood in terms of 
single or multiple scalar fields, the inflaton fields, which 
roll in certain underlying potentials [12]. The action for 
gravity and a minimally coupled scalar field in general 
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relativity is [13] 


represented as 


S^ = - ct^xy/^ 


I - + V{^) 


( 1 ) 


S = — d'^Xy/^ 


R 

'2 




with V{(p) the self-interaction potential of the scalar field. 
We have chosen the natural system of units such that the 
gravitational constant obeys G~^ = Stt, the speed of light 
c = 1 and the reduced Planck constant h~^ = 2y/TT. 

Many different dark energy models were advanced with 
a scalar field responsible for the late-time cosmic acceler¬ 
ation. They include those with a single canonical scalar 
field fp, the quintessence with a non-zero potential [14]. 
Another class is provided by the k—essence models, in 
which the late-time acceleration is driven by the kinetic 
energy term X of the scalar field and the Lagrangian is 
an arbitrary function of (p and X [15]. In its Dirac-Born- 
Infeld subcase [16, 17] L = —V (</>) Vl — 2A, the scalar 
is known as the tachyonic field. When further the poten¬ 
tial y is a constant, the scalar becomes the Chaplygin 
gas, also investigated as a dark energy candidate [18]. 
Some of these models are also unified models for dark 
energy and dark matter, similarly to Refs. [19]. Finally 
we mention the coupled models, where dark energy inter¬ 
acts with dark matter [20]. For a review of the possible 
dark energy candidates see [21]. 

Despite successful scalar field dark energy models, one 
cannot reject a priori the idea that dark energy has a 
more complex field structure. One intriguing possibility 
is to model dark energy as a vector or Yang-Mills type 
field, which may also directly couple to gravity. The sim¬ 
plest action for a Yang-Mills type dark energy model is 
[ 22 ] 


Sv = - j d'^x 
+V{A^) 


where — V^A‘p and is the covari¬ 

ant derivative with respect to the metric, while A^ = 
gU’Aj\a^A'p. In the action given by Eq. (2) dark energy 
is represented by three vector fields. This generalizes an 
Einstein-Maxwell single vector dark energy model, with 
V{A^) representing a self-interaction potential explicitly 
violating gauge invariance. The astrophysical and cosmo¬ 
logical implications of the vector type dark energy models 
have been extensively investigated in Ref. [23]. 

One can construct extended vector held dark energy 
models either, in which the vector held is non-minimally 
coupled to the gravitational held. Such models have been 
proposed, for example, in Ref. [24]. The action for the 
non-minimally coupled vector dark energy model can be 



_ T^a T^afiy 


( 2 ) 


AioA^A^R rjA^A'^Rfju -I- 


( 3 ) 


where fiA is the mass of the massive cosmological vector 
held, and (x*"), pL,v = 0,1,2,3 is its four-potential, 
which couples non-minimally to gravity. By analogy with 
electrodynamics, the dark energy held tensor is dehned 
as = Xfj^Ai, — XvA^. Here w and rj are dimensionless 
coupling parameters. 

Inspired by some condensed matter concepts, a so 
called superconducting type dark energy model was re¬ 
cently proposed [25] . The starting point of this approach 
is represented by the deep connection of the gravitational 
actions given by Eqs. (1) and (2), respectively. Despite 
their different form, they can be interpreted and under¬ 
stood as two limiting cases of a single unihed fundamental 
physical model, describing the spontaneous breaking of 
the U(l) symmetry of the ’’electromagnetic” type dark 
energy. In condensed matter physics such an approach is 
used to describe the superconductive properties of metals 
from a held theoretical point of view [26-28] . 

Hence the basic assumption in the superconducting 
dark energy model is that its ground state is in the form 
of a Bose-Einstein type condensate, with the U(l) sym¬ 
metry spontaneously broken. Consequently, the U(l) 
gauge invariant dark energy appears as a superposition of 
a vector and of a scalar held, the latter corresponding to 
the Goldstone boson. This model is thus equivalent with 
the standard held theoretical description of the supercon¬ 
ducting phase. In [25] the gravitational held equations 
included a non-minimal coupling in the action between 
the cosmological mass current and the superconducting 
dark energy. 

Due to the coupling between dark energy and matter 
current, in the superconducting dark energy model the 
total matter energy-momentum tensor including 

the contributions of both baryonic and dark matter com¬ 
ponent, is not conserved. This non-conservation of the 
matter energy can be interpreted in terms of the ther¬ 
modynamics of irreversible processes as describing mat¬ 
ter creation in an open system. However, in the present 
model particle creation is not the main physical process 
determining the accelerated expansion of the Universe. 

In the framework of a Eriedmann-Robertson-Walker 
homogeneous and isotropic geometry two particular cos¬ 
mological models were considered, corresponding to the 
choice of a pure electric type held, and of a pure mag¬ 
netic type potential, respectively. In both cases the cos¬ 
mological evolution was investigated by numerical and 
analytical methods. As the main result of these studies 
it was shown that in the presence of either the electric or 
the magnetic type superconducting dark energy the Uni¬ 
verse evolves into an exponentially accelerating vacuum 
de Sitter state. 
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Weather this acceleration emerging in the dark energy 
models discussed in Ref. [25] falls in the range given by 
cosmological observations, still needs to be tested. It is 
the goal of the present paper to continue and extend the 
investigation of the superconducting dark energy mod¬ 
els, by confronting them with the most recent supernovae 
type la (SNIa) data set [29], then with the available data 
on the Hubble parameter-redshift relation [30-32]. Such 
tests are nowadays routinely applied to all viable dark en¬ 
ergy models, hence it is a must to confront with them the 
superconducting dark energy models and restrict their al¬ 
lowed parameter spaces accordingly. 

The present paper is organized as follows. In Section II 
we briefly review the mathematical and physical formal¬ 
ism of the superconducting dark energy model. We per¬ 
form the cosmological tests in Section III. We discuss 
and conclude our results in Section IV. 

In what follows we use the Landau-Lifshitz [33] sign 
conventions for the definition of the metric and of the 
geometric quantities. 


II. THE SUPERCONDUCTING DARK ENERGY 
MODEL 


In the present Section we will briefly review the gravi¬ 
tational field equations and the corresponding cosmologi¬ 
cal application of the superconducting dark energy model 

[25] to be used in the sequel. Dark energy is assumed to 
be described by two fields, a vector field (analogous 
to the electromagnetic 4-potential) and a scalar field (j) 
(analogous to the Goldstone boson). Our basic physi¬ 
cal requirement is the breaking of the U(l) symmetry 
of dark energy by an operator of ’’charge” —2e, simi¬ 
larly as for superconductivity (see Section 21.6 of Ref. 

[26] ). Then it follows that the action for a unified scalar- 
vector dark energy theory is invariant under the gauge 
transformations given by A^{x) —>■ A^(x) + d^C{x), and 
4>{x) —>■ 4){x) -I-C'(x), where C{x) is an arbitrary function 
of the coordinates. 

Hence U(l) gauge invariance of the dark energy re¬ 
quires that the fields (<(), A^) should appear only in the 
combinations F^i, = V^Ai, — V^A^ and 

21^ = , (4) 


respectively. 


A. Action and gravitational field equations 


binations F^i, and 21^. The action is: 


S = — dt^Xy/^ 


f + - ^21^ + V (a^) 


a 


+ Lm - 0 ) 


( 5 ) 


where = g^'^a^a^, the notation V (a^) stands for the 
self-interaction potential of the a^ field, A = ^^/2 is re¬ 
lated to the mass of the massive cosmological vector field, 
while Lm {g^i ,, -ip) is the Lagrangian of the total (baryonic 
plus dark) matter, represented by ip. In Eq. (5) = pu^ 

denotes the total cosmological mass current, where p is 
the total (baryonic plus dark) matter density, and is 
the matter four-velocity. We assume that the baryonic 
and dark matter are comoving with the cosmological ex¬ 
pansion. In the action (5) we have also introduced an 
interaction term between the total matter flux and 
the superconducting dark energy gauge invariant poten¬ 
tials A^ — with a representing a coupling constant. 

The action (5) has two important limiting cases. When 
<j) = 0, and A = 0, a = 0, the gravitational action (5) 
gives the single field vector model of dark energy [22, 23], 
a subcase of the action (2). On the other hand, when the 
electromagnetic type potential vanishes, A^ = 0, and 
A = 1, V (a^) = 0, we reobtain the gravitational ac¬ 
tion of the minimally coupled scalar-tensor theory (1), for 
V{(p) = 0. Therefore the gravitational action (5) provides 
a unified framework for the inclusion into the action of a 
dark energy theory of the minimal scalar-vector interac¬ 
tions, under the fundamental assumption of the existence 
of a U(I) type broken symmetry. It is interesting to note 
that the second and third terms in Eq. (5) are similar to 
the corresponding terms in the Stueckelberg Lagrangian 
[34]. Indeed, massive vector fields replace a massless 
vector field and a Goldstone boson in the nonrelativis- 
tic treatment of superconductivity by Anderson [35] and 
the related Englert-Brout-Higgs-Guralnik-Hagen-Kibble 
mechanism of mass generation by spontaneous symmetry 
breaking [36]. 

By varying the action (5) with respect to the metric 
tensor, the scalar field (p and the vector potential A^ we 
obtain the full set of the gravitational field equations for 
the superconducting dark energy model as 




1 

47r 



-k 


A2l^2ti, - ^ - 2 x 

aa2'G(2l')2t^2t,-iu(2l2)5^, , 


( 6 ) 


An action for the superconducting dark energy model 
was proposed in Ref. [25]. We modify that action here 
such that it becomes manifestly gauge-invariant, hence 
the vector and the scalar fields appear only in the com- 


[^212^ (21^) A^] - 

AV^A'^ - = 0, (7) 
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and 


1 

47r 




( 8 ) 


respectively, where we have denoted 


= 


Aa'" + - 2^212^ (21^) 


(9) 


In order to obtain Eqs. (6)-(9) we have used the mathe¬ 
matical identities 


5v 

dE (g^^a^a.) 
5(j) 

dE (5^"a^a,) 
SA^ 


aa2i/ (a^) a^a,<55^^ 

-2d<^2V {^'^) 

2d<^2V g^^’^AJA^, 


and 


-2 J [d^2V (a^) gf^^'A^WM] V^d^x = 

2 Jg^^v, [d^2V (a^) A^] V^d^x, 


respectively, and we have assumed that at infinity the 
variation of the scalar field (j) vanishes. 

In the following we assume that the matter content 
of the Universe consists of a perfect fluid, and therefore 
the energy-momentum tensor of the matter, in¬ 

cluding both the dark matter and the baryonic matter 
components, with energy densities and pressures pdm, 
Pdm and ps and pb respectively, is given by 

^"^'‘T^^ = {p + p)u^Uu-pgfj,^, ( 10 ) 

where p = pdm + Pb and p = pdm + Pb are the total 
thermodynamic energy density and pressure of the mat¬ 
ter components (baryonic and dark). The dark energy 
electromagnetic type tensor automatically satisfies 
the Bianchi identity 

= 0, (II) 

where is the complete antisymmetric unit tensor 

of rank four. 

The conservation equation, including all the compo¬ 
nents of the Universe in the presence of the supercon¬ 
ducting dark energy is obtained as 

+ 

2V^V^ (21^) A'"] -k 29^214 (21^) A“V„Ai. = 0(12) 

In the rest of the paper we restrict our analysis to the 
case of a constant self-interaction potential of the super¬ 
conducting dark energy field, V (21^) = Vb = constant. It 
is interesting to note that in the theory of superconduc¬ 
tors such a potential corresponds to the average energy of 
the ground state of the condensate < eapip^^ip^ > formed 
in the superconductor, with Tp°‘ the wave function of the 
condensate, and Cap the two-dimensional completely an¬ 
tisymmetric unit tensor [26, 27]. 


B. Cosmological models 

For cosmological applications we assume that the met¬ 
ric of the Universe is given by the isotropic and homoge¬ 
neous flat Friedmann-Robertson-Walker metric, 

da^ = dt^ — a^{t) [dx"^ + dy"^ + dz"^) , (13) 

where a{t) is the scale factor describing the expansion of 
the Universe. 

In the following we will discuss two types of supercon¬ 
ducting dark energy models, the electric and magnetic 
type models, respectively. 


1. The electric type (scalar) superconducting dark energy 
model 

In the electric type (scalar) dark energy model we as¬ 
sume that the dark energy vector potential has the simple 
form Afj, = (Ao(t), 0,0, 0), hence only the electric (scalar) 
potential is kept. This dark energy 4-vector potential 
leads to = 0, V/r, g [0,1,2,3], hence both the 
’’electric” and the ’’magnetic” dark energy fields vanish 
and the scalar Ao(t) turns out to be non-dynamical in 
a Lagrangian sense. In other applications of the super¬ 
conducting dark energy models however, notably in the 
static, spherically symmetric case, when Aq = (j){r), a 
pure electric type field is generated. That case is impor¬ 
tant when discussing the Newtonian limit, or the dark 
energy effects on the planetary dynamics and light mo¬ 
tion in the Solar System. Therefore we adopt the unique 
terminology ’’electric type superconducting dark energy 
model” for this case. 

Since = 0, we have = 0 from Fq. (8), then Fq. 
(9) gives 

= (14) 

for a constant potential V (21^) = Vq. Only the com¬ 
bination Aq — <f) occurs in the field equations which can 
be eliminated by —ap) (2A) [25]. Hence the generalized 
Friedmann equations describing the cosmological expan¬ 
sion in the electric type superconducting dark energy 
model become 

3^2 = p - ^p^ + Vq = p + pde , (15) 

2H + 3i7^ = -p+ —p^ + Vq = -p- Pde , (16) 

where we have denoted 

q ,2 

Pde = + ^o,Pde = ~ , ( 17 ) 

and F[ = a/a is the Hubble parameter. From Fqs. (15 ), 
then Fq. (16), we have 
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Then from the derivative of Eq. (15) with respect to time 
and eliminating H by Eq. (18), we find 

■ Ix'’) + P “ 

For dust {p = 0) this reduces to 

(p + 3ifp)(^l-|^p) =0, (20) 

giving a non-constant solution p oc a~^ expressed in 
terms of redshift {z = ao/a — 1): 

p{z) = po {I + zf , (21) 

where the subscript 0 denotes the present value. Finally 
the Hubble parameter is given by 

rp' 

3H^{z) = po(l + - 

Despite of a ^ 0, the four-divergence of the dust’s 
energy-momentum tensor vanishes in Eq. (12), and the 
matter density satisfies a continuity equation. However a 
non-vanishing coupling between the matter and the fields 
A and (j) leads to the excitations of the dark energy as 
can be seen from Eq. (14). 


2. Magnetic (vector type) superconducting dark energy 
models 


In the following we call magnetic type supercon¬ 
ducting dark energy models the cosmological models 
in which the 4-vector potential is given by = 
{0,Ai{t),A2{t),A3{t)), hence by the magnetic (3-vector) 
potential. For an isotropic expansion the components of 
Afj^ satisfy the condition Ai{t) = A 2 {t) = A^iJ,) = A(t). 
However, this model has the interesting property that 
because of the choice of the potential and of the space- 
time metric, and its symmetry properties, the ’’magnetic 
field” identically vanishes^ but the electric field is non¬ 
zero. Again we note that in the static, spherically sym¬ 
metric case a proper dark energy magnetic field emerges 
from the vector potential A = A{r). Therefore, for the 
uniqueness of the terminology, we named this case ac¬ 
cording to the magnetic type potential as magnetic type 
(vector) superconducting dark energy model. 

With the constant value V (21^) = Vo and homoge¬ 
neous, isotropic 4-vector potential, the field equation (8) 
gives J* = 0 which due to Eq. (9) leads to 




ap 

^ ■ 


(23) 


By this expression cj) can be eliminated from the remain¬ 
ing equations. The generalized Friedmann and the en¬ 
ergy conservation equations take the form 




3 A^ 
Stt aP 




4^0 , 


(24) 


2H + ?,H^ = -p+—p^ 


1 

o-2" ■;r —T + Eq 

Stt 2 


(25) 


A-I-HA-I- 47rAA = 0 . 


From Eqs. (24) and (25) we have 


2H = -p 



1 42 

27r aP 


(26) 


(27) 


Then from the time-derivative of (24) and eliminating H 
and A by Eqs. (27) and (26), respectively, we find the 
same evolution equation (19) for the matter density than 
in the electric-type dark energy model case. 

The dust satisfies a matter conservation separately, 
however Eq. (23) tells that the matter is a source of 
the scalar field cj) for a non-vanishing value of a . 


III. TESTING THE SUPERCONDUCTING 
DARK ENERGY MODEL WITH SUPERNOVAE 
AND HUBBLE PARAMETER DATA SETS 


The confrontation of a cosmological model with the 
Union 2.1 SNIa data set [29] through a y^-test, as de¬ 
scribed in Refs. [17] relies on fitting with the computed 
luminosity distance (dz, )-redshift relation. For a flat 
Friedmann universe a dimensionless luminosity distance 
dh = HodL is given by 


d dL 1 

dzl -\- z ff ’ 


(28) 


where z is the redshift, H (z) is the Hubble parameter, 
Hq = H {z = 0) and H {z) = H (z) /Ho- Hereafter the 
subscript 0 always denotes the present value of the cor¬ 
responding quantity. 

Using the data set [30-32] on the Hubble parameter- 
redshift relation we compute 


2 ' \Hth (Zi) HqIjs izi)^ 

xh = 2^ ^2 


(29) 


where Hth (zi) and idobs (zi) are the Hubble parame¬ 
ter values at redshifts Zi, predicted by the cosmological 
model, and determined from the observations, respec¬ 
tively, while CTi is the scattering in Hobs (zi). We note 
that a subset of the data set on the Hubble parameter- 
redshift relation was used to emphasize a tension with 
the ACDM model (by computing the two-point Omh^ 
function) [37]. In the following we perform a test of 
the superconducting dark energy model with the com¬ 
bined SNIa and Hubble parameter data set by calculat¬ 
ing xlNia+H = xlNia+XH, where xIniu is the x^-value 
from the confrontation with the SNIa data set. 

In our analysis we will investigate only the late cos¬ 
mological evolution. Then the mixture of the baryonic 
and the dark matter component has negligible pressure, 
so that in the field equations we can take p = 0. 
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SNIa test 




FIG. 1: (Color online) The test of the electric dark energy 
model with the supernovae data set is shown. The color scale 
represents the values on the parameter space 
Within the gray area the condition (34) is violated. The light- 
blue and the green lines bounds the Icr and 2cr confidence 
regions, respectively. 


A. Cosmological tests of the electric dark energy 
model 


In this model the dark energy vector potential is 
parallel with the four velocity of the comoving mat¬ 
ter. The evolution of the Hubble parameter follows from 
Eqs. (21) and (22), and is given by 


{z) — flm (1 + ■ 2 ) + i^DE ~ Hy) (1 z) fly , 

(30) 

where flm and floE are the present day density param¬ 
eters of the dust and dark energy, respectively: 


H777, — 




z=0 


(31) 


SNIa+H test 



a 

FIG. 2: (Color online) The test of the electric dark energy 
model with the combined supernovae-l-Hubble parameter data 
set is shown. The color scale represents the x^ values on the 
parameter spaces {Qv,^m) (upper panel) and (5,flm) (lower 
panel). Within the gray area the condition (34) is violated. 
The light-blue and the green lines bounds the Icr and 2a con¬ 
fidence regions, respectively. 


^DE = 1 — , (32) 


and 


a = 


8A 


Hy = 


^0 

3iJ2 


(33) 


are dimensionless constants. We constrain the model pa¬ 
rameters by imposing that the universe continuously ex¬ 
panded after the Big Bang Nucleosynthesis (BBN). Since 
{z) > 0, we have the following constraint on the pa¬ 
rameter space, 


(I -b zbbn)^ 

(1 -b ZbBn)^ 


(34) 


The value of 5cr is of the order of 10“^^ for reason¬ 
able cosmological parameters, thus we investigate only 
the negative a range. 


In the test of electric dark energy model the Hubble 
constant is fixed by Hq = 67.74 km/sec/Mpc in agree¬ 
ment with the measurements of Planck spacecraft on the 
anisotropies of the cosmic microwave background radia¬ 
tion (CMB) (see the last column of Table 4 in [38]). The 
results of the tests with the Union 2.1 SNIa data set are 
shown on Fig. 1 in the parameter space (Bm,By). The 
color scale represents the values and the light-blue and 
the green lines bounds the Icr and 2cr confidence regions, 
respectively. The grey area shows the region where a 
bounce occurs between the redshift values 0 and Zbbn, 
i.e. where a > ctcr ■ Good fit with SNIa data set is 
obtained along a narrow inclined stripe, which neverthe¬ 
less contains the ACDM limit (5 = 0), with the widely 
accepted values of ~ 0.3 for the dark plus baryonic 
matter sector and By Ri 0.7 for dark energy. 

The test with the combined SNIa and Hubble param¬ 
eter data set further narrows the Icr and 2cr confidence 
regions, as shown on the upper panel of Fig. 2. On the 
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SNIa+H test 


a 



FIG. 3: (Color online) The test of the magnetic dark en¬ 
ergy model with the SNIa (upper panel) and the combined 
SNIa+Hnbble parameter (lower panel) data sets are shown. 
The color scale represents the values on the parameter 
spaces (A,5). The light-bine and the green lines bounds the Icr 
and 2a confidence regions, respectively. The fixed parameters 
of the model are flm = 0.3088, Ao = 0.05, wde\z=o = —0.99, 
Ho = 67.74 km/sec/Mpc and dA/dz\z^o > 0. The gray area 

^dA/dz^ 


gives the region where 


^0 would be negative for 


the chosen parameter values. 


lower panel the result of the test with the combined data 
set is also represented in the parameter space (q',^^)- 
The full fit in the 3-parameter space {a,Q.rmHo) in 
the case of electric type model gave local minima, among 
which the LCDM limit was the best fit (the extra pa¬ 
rameter 5 = 0). Thus the electric type models favour 
the ACDM models. 


B. Cosmological test of a magnetic dark energy 
model 

In this model the dark energy vector potential is per¬ 
pendicular to the four velocity of the comoving matter. 
In accordance with the isotropy of the assumed space¬ 


time geometry, all spatial components of the vector po¬ 
tential are the same A (z). From Eqs. (24) and (27), the 
Hubble parameter and its derivative are given by 


(z) 


l-(l + zf 



(1 + ^) + 


-5H^ (1 + zf + X{l + z f , (35) 


H (z) dH (z) 
1 + z dz 


-Xlyri (1 -f z) 


35fl^ (1 + zf + AA^ 


+2(l + zfff^ (z) 



(36) 


where A{z) = A (z) /v^S^oo and 

47rA 


A = 


772 


(37) 


The evolution of the dark energy vector potential is gov¬ 
erned by Eq. (26), and can be formulated as a function 
of the redshift: 


d^A 


dA ~ A 

K ——h A- 2 —;— 

dz (l + z) 


= 0 , 


where 


(38) 


K = 3»m(l + z) 

2 H2 JJ2 


+2{l + zf 


dA 

dz 


+ A (1 -l- z) 


Af 


(39) 


Since Ho = 


(^A/dz^ 


\z=0, 


1, Eq. (35) gives a constraint among 
Ao, rtm, Aly, 5 and A. Taking into ac¬ 


count the Hubble constant, the evolution of the super¬ 
conducting dark energy dominated Universe depends on 
six parameters. Since Eq. (35) is a first integral of the 
system (36) and (38), it is used to check the accuracy of 
the numerical solution. 

The Hubble constant, and the cosmological param¬ 
eter measuring the sum of baryonic and dark matter 
components, the matter density parameter are fixed as 
Hq = 67.74 km/sec/Mpc and = 0.3088, respectively. 
These values emerge from the confrontation of the six- 
parameter ACDM model with the observations of the 
Planck satellite (last column of Table 4 in [38]). The 
parameter wde, giving the equation of state of dark en¬ 
ergy at present time, is chosen as 


wde\,,=o 


i {dA/dz'^ 1^=0 + 

aW^-[dA/dz^ [2=0 - AAq - Dy 
-0.99, 


( 40 ) 
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FIG. 4: (Color online) The same as on Fig. 3 but the fixed 
parameters are given by flm = 0.0486, Aq = 0.05, wde\z=o ~ 
—0.855, Ho = 70 km/sec/Mpc and dA/dz\z^o > 0. 


which is consistent with the Planck satellite’s measure¬ 
ments [38]. We fix Aq as 0.05, and choose dA/dz\z=o > 0, 
such that the function A {z) is decreasing at present. 


Then the initial value 


as 


(^dA/dz^ 


can be expressed 


z=0 


^ |z=o = 3{1 — wde)cA}^ 

+3 (1 -l- wde) G -I- (1 -I- Swde) -^^0 I (41) 

where 

G = 1 — ~ •^^0 ■ (42) 

We note that Eq. (41) admits negative values for 
for certain values of Aq , Qm, wde\z= 0 j 

2=0 

5, and A, which need to be excluded. 

The results of the test of the magnetic dark en¬ 
ergy model with the SNIa data and with the combined 




SNIa-l-Hubble parameter data sets are shown in Fig. 3 
in the parameter space ^A,5^). We note that a bounce 

would occur between the redshifts 0 and zbbn for a > 0 
(except with the infinitesimal positive values as in the 
electric dark energy model), and therefore the Figure 
shows only the allowed parameter range 5 < 0. 

We may also ask if the magnetic type superconduct¬ 
ing dark energy fluid can mimic the effects of both dark 
energy and dark matter, at least regarding the data set 
on SNIa and Hubble parameter? In other words, is it 
possible to have a unified dark energy and dark matter 
model? To answer this question we fix the parameter 
to the value 0.0486 of the baryonic matter alone [38]. In 
order to reach good fit with the data, other cosmologi¬ 
cal parameters could also change. We show the result of 
the test on Fig. 4 for wde\z=o = —0.855 and Hq = 70 
km/sec/Mpc as an example. As in the previous case, a 
bounce between the redshifts 0 and Zbbn is avoided for 
5 < 0. It is interesting to observe that the la and 2a 
confidence regions are larger in this case. 

The values of the free parameters of the model which 
are preferred by the cosmological observations include 
a low value of the coupling between the matter current 
and the fields (a ~ 0), while the mass of the dark energy 
vector field increases from A = 0.01663 Hq to A = 0.04584 
Hq when it is also supposed to account for dark matter. 


1. Best fit parameters 


TABLE I: Some local minima of x^. The coupling constant 
5 = 0 and hioo is defined by Ho = lOOhioo km/s/Mpc. 


Xmin 

Wde 

A 

Aq 

hioo 

578.76 

0.148 - 0.983 4.568 

- 0.030 

0.698 

578.99 

0.246 - 0.916 

0.175 

0.528 

0.697 

579.57 

0.268 - 0.997 

- 0.009 

- 0.373 

0.696 


In this subsection we fix none of the six parameters of 
the model. It is well-known that determining the global 
minimum of a function of six variables proves tremen¬ 
dously difficult due to the fact that any known algorithm 
for determining such a minimum could easily run into a 
local minimum. One of the most widely applied heuris¬ 
tic methods for a large parameter space is the Simulated 
Annealing Method. We have applied this procedure as 
implemented in Mathematica [39] in order to find the 
minimum of the chi-squared function. The Simulated 
Annealing method lead to a multitude of local minima, 
suggesting severe degeneration in the parameter space, 
or in other words, elongated 1-sigma contours in the full 
parameter space. The obtained smallest local minimum 
is = 578.76. A x^-difference test, as described in Ref. 
[40], which showed that the fitting to the data sets is of 
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FIG. 5: (Color online) The evolutions of H (z) in the magnetic 
type dark energy model (green, blue, red) and in the ACDM 
model (black). Each curve represents local minima for a-t 
points listed in Table I. On the upper panel the data set for 
Hubble parameter is also shown. 


closely similar quality for the investigated models as com¬ 
pared to the ACDM model. In other words, the ACDM 
model, arising as a special case in the family of models 
investigated, belongs to the Ict region of the parameter 
space. 

Among the local minima we have found evolutions al¬ 
most indistinguishable from ACDM, but also very differ¬ 
ent evolutions, with comparable fit. We represented the 
following local minima in Table I and on Figs. 5 and 6: 
i) the evolution belonging to the smallest local minimum 
found (green curve), ii) the closest evolution to ACDM 
(red curve), and in) an evolution in between the two 
(blue curve). The evolutions of the Hubble parameter 
in the magnetic model (green, blue and red curves) and 
for the ACDM model (black curve) are shown on Fig. 



FIG. 6: (Color online) The evolutions of A {z) related to the 
curves shown on Fig. 5. 

5. On the upper panel the evolutions are represented for 
low redshifts and the data set for the Hubble parameter 
are also given. The lower panel shows the evolutions for 
high redshifts. The evolutions of A (z) are represented on 
Fig. 6. In the distant past A (z) approximates a constant 
value, which is a typical behaviour in this model. 


IV. CONCLUDING REMARKS 

In the present paper we have considered the cosmolog¬ 
ical tests of an electromagnetic type dark energy model 
introduced in Ref. [25] , in which the electromagnetic 
gauge invariance of the dark energy is spontaneously bro¬ 
ken. From a physical point of view the superconducting 
dark energy model considered in our analysis is a two- 
field model, and extends general relativity into a scalar- 
vector-tensor gravitational theory. Alternatively, it can 
be seen as a unified scalar - vector field dark energy 
model, in which the scalar field (j) and the vector field 
appear in the gauge invariant, massive combination 
— V^(t>. From the point of view of cosmological 
applications, the present model may be used for realiza¬ 
tions of early time inflation, with a possible explanation 
of the large scale CMB anomalies, or, as done in the 
present paper, for describing the late time acceleration 
of the Universe, and as a specific model of dark energy. 

An important theoretical test of gravitational models 
is the analysis of the cosmological perturbations. The 
stability results also have important observational im¬ 
plications. In vector-tensor models with action given by 
Eq. (2), or in more general models with non-minimal cou¬ 
pling between vector fields and gravity, the analysis of the 
stability proved to be a challenging task. In contrast to 
the scalar field gravitational models, the linearized the¬ 
ory of vector-tensor models contains couplings between 
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scalar, vector and tensor modes [41]. It turns out that 
the perturbations decouple only in the ultraviolet limit 
[41], which, on the other hand, opens the possibility of 
the quantization of the models. Note that the mixing of 
the different modes in vector - tensor models may lead 
to observable correlations of the scalar and tensor modes 
in the CMB. 

In [42] it was shown that a number of vector-tensor 
gravitational models contain instabilities in the form of 
ghosts, or unstable growth of the linearized perturba¬ 
tions. These instabilities can be related to the longitu¬ 
dinal vector polarization modes presents in the theoret¬ 
ical models. The presence of ghosts or tachyons in the 
early time/small wavelength regime of the vector-tensor 
cosmological models are an indication that the vacuum 
of the model is unstable. The instabilities found in [42] 
growth up to horizon crossing, thus indicating that the ef¬ 
fective cosmological constant A is comparable, or smaller 
than the Hubble scale, implying a different cosmological 
phenomenology. An analysis of perturbations in vector 
field models, with a special emphasis on inflation, was 
performed in [43], where the long wavelength limit of 
perturbations in small fields inflation and their linear sta¬ 
bility was discussed. It turns out that the gravitational 
waves are unstable in large fields models. However, as a 
main result of this investigation, it was shown that one 
can use the small fields approximation in vector-tensor 
models as an effective theory for describing cosmological 
evolution. 

The Hamiltonian stability and the hyperbolicity of vec¬ 
tor field models involving both a general function of the 
Faraday tensor and its dual, as well as a Proca type po¬ 
tential for the vector field, were considered in [44]. The 
main result of this investigation is that the considered 
vector models do not satisfy the hyperbolicity conditions, 
which require that the Cauchy problem is well posed. 
However, as shown in [45], in vector-tensor field theories 
with non-canonical kinetic terms, violations of hyperbol¬ 
icity may not be present around the cosmological back¬ 
grounds. Therefore vector models involving new cosmo¬ 
logical dynamics and no causal pathologies can be found. 

The problems of the stability of the linear perturba¬ 
tions and of the hiperbolicity is also of fundamental im¬ 
portance for the superconducting dark energy model. A 
full analysis can be done by considering the perturbations 
of the field equations Eqs. (6), (7) and (8), respectively, 
and analyzing their stability. However, we would like to 
point out that the superconducting dark energy model is 
explicitly constructed as a stable ground state of the U(l) 
gauge invariant scalar-vector field configuration. Hence 
we expect that this configuration is stable at least with 


respect to small linear perturbations, as is also the case 
for other similar vector-tensor models [45]. 

We have considered the observational tests for two dis¬ 
tinct classes of models, whose physical and cosmological 
properties are determined by the choice of the dark mat¬ 
ter electromagnetic potential A^. The first model cor¬ 
responds to an electric type dark energy potential, with 
Aq (t) the only non-zero component, while in the second 
case the space-like components of are different from 
zero, with A^ = {0, A{t), A{t), A{t)). For simplicity in 
both cases we assumed a constant self-interaction poten¬ 
tial for the electromagnetic type dark energy with U(l) 
broken symmetry. 

In the case of the electric model good fit with SNIa 
data set is obtained along a narrow inclined stripe in the 
parameter space (0^, Hy), which nevertheless contains 
the ACDM limit (5 = 0), with the widely accepted values 
of Qm ~ 0.3 for the dark plus baryonic matter sector 
and fly Ri 0.7 for dark energy. The other points on 
this admissible region represent essential deviations from 
the ACDM model (the black region on the lower panel 
of Fig. 2), but they align along a narrow stripe (upper 
panel of Fig. 2) with similar slope as the one given by 
the SNIa test of the ACDM model, indicating that the 
dominant part of the superconducting dark energy is still 
represented by a cosmological constant type contribution, 
supplemented by a time-evolving component. 

In the magnetic case the cosmological test select either 
i) parameter ranges of the superconducting dark energy 
allowing for the standard baryonic plus dark matter con¬ 
tent flm ~ 0.3 or ii) a unified superconducting dark mat¬ 
ter and dark energy model, including additionally only 
the baryonic sector ~ 0.05. A glance on Figs. 3 
and 4, giving the admissible parameter spaces in the two 
cases shows that 

a) the mass term A in the action has to increase in 
order to incorporate dark matter in a unified model, and 

b) the cosmological data is best matched for a « 0, 
which decouples the matter current from both the scalar 
and vector sectors of dark energy. 
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